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Abstract 

We study the impact parameter dependence of inelasticity in the framework of an updated geometrical model for multiplicity 
distribution. A formula in which the inelasticity is related to the eikonal is obtained. This framework permits a calculation of the 
multiplicity distributions as well as the inelasticity once the eikonal function is given. Adopting a QCD inspired parametrization 
for the eikonal, in which the gluon-gluon contribution dominates at high energy and determines the asymptotic behavior of 
the cross sections, we find that the inelasticity decreases as collision energy is increased. Our results predict the KNO scaling 
violation observed at LHC energies by CMS Collaboration. 
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1 Introduction 

Inelasticity is defined as the fraction of the available en- 
ergy relesead for multiple particle production in inelastic 
hadronic interactions. The remaining part of the incident 
energy is carried away by the participant's remnants, so 
called leading particles. The energy dependence of in- 
elasticity is a problem of great interest from both the- 
oretical and experimental standpoints [1]. However, the 
experimental data are scarce and, on the theoretical side, 
the existing models are largely in conflict with each other 
even in explaining a simple aspect as the center-of-mass 
energy dependence of the inelasticity [2] , [3] . For exam- 
ple, the decrease in inelasticity with energy is advocated 
by some authors, while others believe that the inelas- 
ticity is an increasing function of energy [1]. Hence the 
problem remains unsolved. Naturally, multiplicity dis- 
tributions are connected to inelasticity ones, so one can 
study multiplicity distribution features in order to de- 
rive information on the inelasticity behavior. Following 
this way, we have updated an existing phenomenologi- 
cal procedure [4], referred as Simple One String Model, 
which allows simultaneous description of several experi- 
mental data from elastic and inelastic channels through 
the Unitarity Equation [4]. Thus, based on the Sim- 
ple One String Model formalism, able to describe the 
charged multiplicity distributions from ISR to Collider 
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energies (30.4 - 900 GeV), we have computed the impact 
parameter dependence of inelasticity at fixed center-of- 
mass energies, y^. We have also inferred information on 
energy dependence of inelasticity. The plan of the pa- 
per is as follows. In the section 2 we present the basic 
formalism of the Simple One String Model and the pre- 
dictions for pp/pp overall multiplicity distributions com- 
pared with the experimental data. In Section 3 we apply 
the theoretical framework computing the inelasticity of 
hadronic reactions. The final remarks are the content of 
Section 4. 



2 Simple One String Model for Multiplicity Dis- 
tributions 

The Simple One String Model has been discussed in ref- 
erences [4], [5] and in order to define the notation and 
also update the model, we shall review the main points 
here. We work in impact parameter space, 6, and to 
guarantee unitarity, the inelastic cross sections, aim is 
calculated via the relation: 

a^n{s) = J dHGin{s,b), (1) 

where 

Ginis,b)^l~e-^''''^''''^ (2) 
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is the Inelastic Overlap Fimction. In this work we have 
update the model adopting the complex eikonal func- 
tion xlpis,b) = XR{s,b) + ixi{s,b), from Ref. [6], as 
will be discussed in subsection 2.1. In Ref. [4] has been 
adopted the Henzi Valin parametrization for Gin{s,b). 
The probabilities of n particle production, namely mul- 
tiplicity distribution, P„, is the most general feature of 
the multiparticle production processes [7] and measure- 
ments of charged particle multiplicity distributions have 
revealed intrinsic features in pp / pp interactions [8] . The 
multiplicity distribution is defined by the formula [9] 



Puis) = 



cr„(s) 



T,n=0 (^n{s) (Tin{s)' 



(3) 



where cr„ is the cross section of an n-particle process (the 
so-called topological cross section). The charged mul- 
tiplicity distribution, in the impact parameter formal- 
ism, may be constructed by summing contributions com- 
ing from hadron-hadron collisions taking place at fixed 
impact parameter. In this way, the idea of a normal- 
ized multiplicity distribution at each impact parameter 
b is introduced [10]. Thus the multiplicity distribution is 
written as 



Puis) 



a„(s) !d^bG,n{s,b) 



Jd%G,n{s,b) 



(4) 



where the topological cross section (j„ is decomposed 
into contributions from each impact parameter b with 
weight Gin{s, b). In the original formulation [4] the quan- 
tity em brackets scales in KNO sense, and the Eq. (4) 
can be written as 



Puis) 



[Ml, GUs,b) 
J " " <n(s,b)> 



jd?bGUs,b) 



(5) 



where < n(s, 6) > is the average number of particles 
produced at b and \fs due to the interactions among 
hadronic constituents involved in the collision and, in 
this model, < n(s, b) > factorizcs as [4] 



<n(s,6) >=<iV(s)>/(s,6), 



(6) 



where < N{s) > is the average multiplicity at and 
f{s,b) is the so called multiplicity function. Similarly to 
KNO, is introduced, for each b, the elementary multi- 
plicity distribution 



(1) 



=< n(s, 6) > 



CT„(s,&) 



<n{s,b)>J ' crin{s,b)' 

Thus, with Eqs. (6) and (7), Eq. (5) becomes 



(7) 



f G,„is,b) , (1) / z \ 
Jd%Gin{s,b) 



(8) 



where <l)(s, z) N{s) > P„(s) and z = n/ < N{s) >. 
Here z represents the usual KNO scaling variable. Now, 
to obtain the multiplicity function /(s, 6) in terms of the 
imaginary eikonal X/, it has been assumed that 

(1) the fractional energy that is deposited for par- 
ticle production in a collision at b is proportional 
to XI- 

V7 = 0{s)xi{s,b). (9) 
The physical motivation of this equation is that 
eikonal may be interpreted as an overlap, on the im- 
pact parameter plane, of two colliding matter dis- 
tributions [10]; 

(2) the average number of produced particles depends 
on the energy in the same way as in e^e+ an- 
nihilations, which is approximately represented by 
a power law in ^/s [4] 



< n(s, 6) >= 7 — 

V Sn 



(10) 



where So=l GeV^. In Ref. [11] a power law energy de- 
pendence of multiplicity in both pp and Pb + Pb colli- 
sions has been analyzed based on the gluon saturation 
scenario. Now, combining Eqs. (9), (10) and (6), we have 



f{s,b) 



< N{s) > 



2A 



{Xi{s,b)f^=^{s)[xi{s,b)f^, 
(11) 



with ^(s) determined by the usual normalization condi- 
tions on $ [4] and that serves to determine ^(s) as an 
energy dependent quantity, explicitly [4] 



J dHG,nis,b) 



Jd^bGUs,b)[xi{s,b)] 



2A • 



(12) 



Thus, adopting an appropriate parametrization for Gj„ 
and tp^^\ as well as an adequate value for A, we can 
test the formalism embodied in Eq. (8) and (12), making 
direct comparisons with multiplicity distribution data. 
In the following, we will discuss the results obtained in 
the context of our updated model. 

2.1 Inputs and Results on Multiplicity Distributions 

The Simple One String Model is based on the idea of 
multiparticle creation due to the interactions between 
hadronic constituents in collisions taking place at b. It is 
assumed that in parton-parton collision there is forma- 
tion of a string, in which probably one qq has triggered 
the multitude of the final particles. In previous analysis 
the authors [4] considered the experimental data avail- 
able on e~e~^ annihilations as possible source of informa- 
tion concerning elementary hadronic interactions and, in 
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this work, wc borrow their rcsuhs. Specifically, assuming 
a gamma distribution normalized to 2, 



i^^'\z) = 2^^z''-^e-''^, (13) 

experimental data on e~e+ multiplicity distributions 
were fitted, obtaining fc=10.775 ± 0.064 {x^/Ndf = 
2.61). Also, the average multiplicity data in e^e+ anni- 
hilations were fitted by Eq. (10), giving .4=0.258 ± 0.001 
and 7=2.09 {x^/Ndf = 8.89) in the interval 5.1 GeV 
< ^< 183 GeV, and ^=0.198 ±0.004 {x^/Ndf = 1-7) 
for the set in the interval 10 GeV < < 183 GeV, 
respectively [4]. In the analysis done in Ref. [11] the 
value of A=0.11 was obtained within the gluon satu- 
ration picture. The difference between the values of A 
is, probably, associated with the different sets of ex- 
perimental data used in each analysis. We recall that 
the value of ^=0.11 was obtained by using experimen- 
tal data for average multiplicity of hadrons in the gluon 
and quark jets in e~e+ annihilation, in the interval of 
the jet energy between 0.6 ^ 32 GeV [11]. At the end 
of the section the One String Model formalism will be 
tested using the three A values, above mentioned. Now 
is needed a parametrization for the eikonal function, 
and we have adopted the QCD- inspired complex eikonal 
from the work of Block et al. [6] in which the eikonal 
function is written as a combination of an even and odd 
eikonal terms related by crossing symmetry ~)^{s, b) = 
X^{s,b) ± X (s,&). The even eikonal is written as the 
sum of gluon-gluon, quark-gluon and quark-quark con- 
tributions: 

X+(S, b) = Xqq{s, b) + Xqg{s, 6) + Xgg{s, b), (14) 

while the odd eikonal, that accounts for the difference 
between pp and pp, is parametrized as 

X-(s,6) = C-^^e-/4W^(6;M-). (15) 

The various parameters and functions involved in last 
two expressions are discussed in Ref. [6]. By fixing the 
value of A=0.258, V'^^^ given by Eq. (13) with A;=10.775, 
adopting Gm from analysis by Block et al. and observing 
that ^(s) is obtained by Eq. (12), we have computed the 
overall multiplicity distributions arising from pp/pp col- 
lisions at energies' 52.6, 200, 546 and 900 GeV. The the- 
oretical curves are shown in Figs. 1, 2, 3 and 4 together 
with the experimental data. The curves at pp-ISR 52.6 
GeV and CERN pp Collider 546 GeV shows excellent 
agreement with the data. Figs. 1 and 3 respectively. At 
energies Vs=200 and 900 GeV, also in CERN pp Col- 
lider, the agreement with data seems reasonable since the 
curves agree with experimental points for z' > 1, (high 
multiplicities). Figs. 2 and 4 respectively. In view of the 
recent results on multiplicity distributions, in pp colli- 
sions at •\/s=0.9, 2.36 and 7 TeV at the Large Hadron 
Collider (LHC), reported by CMS Collaboration [15], 
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Fig. 1. Overall scaled multiplicity distribution data for pp 
at ISR energy [12], compared to theoretical prediction using 
the Simple One String Model, Eqs. (8) and (12). 
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Fig. 2. Overall scaled multiplicity distribution data for pp 
at Collider energy [13], compared to theoretical prediction 
using the Simple One String Model, Eqs. (8) and (12). 

we have also computed the overall multiplicity distribu- 
tions that the One String Model predicts at LHC ener- 
gies. The results are shown in Fig. 5 and we can see vio- 
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lation of KNO scaling, in qualitative agreement with the 
result obtained by CMS Collaboration in pseudorapid- 
ity interval of \ri\ <2.4, as discussed in Ref. [15]. As men- 
tioned, two fit values of A=0.258 and A=0.198 have been 
obtained in the previous study [4] , the first one giving a 
better account of lower energy data whereas the second 
one higher energy data. As before, we have computed the 
corresponding hadronic multiplicity distribution by fix- 
ing both the gamma parametrization for Eq. (13), 
and the complex eikonal, Eqs. (14) and (15), and con- 
sidered the two parametrizations for the average multi- 
plicity, ~ sO-258 (^(s)=1.424) and~ s° i98 (^(s)=1.348). 
The results at 546 GeV are shown in Fig. 6 and, for 
A=0.198, we can see the disagreement of the theoreti- 
cal curve when compared to the data. As pointed out 
in [4] the parametrization ^ gO i^s j^j-j^gg information 
from data at high energies, while the parametrization 
~ jg jjj agreement with data at smaller energies. 

However, the information from the e~e+ average multi- 
plicities at high energies does not reproduce the overall 
multiplicity distribution. Hence, by using A=0.256 the 
output seems to be more consistent with data. In addi- 
tion, there is no evidence of gluon saturation at CERN 
pp Collider 546 GeV, however and as a pedagogical ex- 
ercise, we have also computed $ considering the value 
of A=0.11 (C(s)=1.211), Fig. 6. We would like empha- 
size that when the formalism is applied, considering the 
three values of ^ (0.258, 0.198 and 0.11) at energies 52.6, 
200 and 900 GeV, the results are essencially the same 
obtained in Fig. 6. We have expressed $ in terms of mod- 
ified the scaling variable z' = n — No/ <n — No > with 
iVo=0.9 representing the average number of leading par- 
ticles [4]. 



3 Inelasticity 

The concept of inelasticity is essential since it defines 
the energy available; for particle; prochiction in high en- 
ergy hadronic and nuclear collisions. However, the im- 
pact parameter dependence of the inelasticity is a prob- 
lem unsolved. In theoretical works, it is quite natural 
to assume that the multiplicity distribution and inelas- 
ticity are connected. Indeed, some authors has defined 
multiplicity distributions in terms of inelasticity as [1], 
fl6l 
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Fig. 3. Overall scaled multiplicity distribution data for pp 
at Collider energy [14], compared to theoretical prediction 
using the Simple One String Model, Eqs. (8) and (12). 



Puis) 



P{n\K)P[K{s)]dK, 



(16) 




where P[K{s)] is the inelasticity distribution and 
P{n\K) is the probability of the production of n par- 
ticles at the given inelasticity K. Thus, based on the 
connection between multiplicity and inelasticity, we 
have explored the parametrization of the One String 
Model formalism and computed the impact parameter 
dependence of inelasticity, as will discuss in next two 
subsections. 



Fig. 4. Overall scaled multiplicity distribution data for pp at 
Colider energy [13], compared to theoretical prediction using 
the Simple One String Model, Eqs. (8) and (12). 
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observed [16], [17], [18]. For example, the average mul- 
tiplicities in pp/pp and e~^e~ collisions become similar 
when comparisons are made at the same effective energy 
for hadron production. In pp/pp collisions the effective 
energy for particle production, E^ff, is the energy left 
behind by the two leading protons 

1 + Eleading,2), (17) 

or 

^eff = {Vs)pp — '2Eieading, (18) 

in the case of symmetric events. (We recall that (v^)pp 
and represents both the center-of-mass energy, how- 
ever, in this subsection the notation {\/s)pp is helpful to 
differentiate that from (v^)e+e- )• In e+e" collisions the 
effective energy for hadron production coincides with to- 
tal center-of-mass energy of the beam 



E. 



eff 



2Ei 



beam' 



(19) 



Thus, the same equivalent energy for both pp/pp and 
e'^e^ collisions can be written as 



Fig. 5. Theoretical predictions for overall multiplicity distri- 
bution by using the Simple One String Model, Eqs. (8) and 
(12), at LHC energies. 



{Vs)pp — 2Eleading = Egff = (v^)e+e-- 



(20) 



For the quantitative estimation of the inelasticity, K, we 
can use the definition [1] , [2] 
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In the following, we will explore the One String Model 
formalism, able to describe the multiplicity distributions 
in wide interval of energy (30.4 - 900 GeV), to obtain 
information about inelasticity. 

3.2 Computation of Inelasticity 

The Eq. (9) is a key point of the formalism. Physically, 
it corresponds to the energy for hadron production de- 
posited at b, due to the interactions among hadronic 
constituents involved in the collision. Thus, and as dis- 
cussed in last section, the fractional energy, \/s' (Eq. 
(9)), and the effective energy for hadron production, 
Epff (Eq. (21)), represents both the same physical quan- 
tity {W = Eeff). Now, by using the Eq. (9), let us write 
the inelasticity, Eq. (21), as a function of ^/s and b as 



Fig. 6. Overall scaled multiplicity distribution data for pp at 

Collider energy [14], compared to theoretical prediction using 
the Simple One String Model, Eqs. (8) and (12), considering 
three different values of A, Eq. (10). 

3.1 Similarities between pp/pp and e'^e~ Collisions 

The idea of a universal hadronization mechanism is not 
new and similarities between both processes were indeed 



2K{s,b) 



(^A), 



(22) 



The factor 2, in the above equation, is due to the fact 
that the multiplicity distributions data are normalized 
to 2. However, we can not calculate the K[s,b) until the 
value of /3(s) is known. To estimate /3(s) we note that 
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the parameter £,{s), which is introduced in Eq. (11), is 
related with /3(s) by 
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< N{s) > 



n 2A 



(23) 



We recall that < N{s) > is the average multiplicity at 
^/s. By using the values of ^ = 0.258, 7=2.09, discussed 
in subsection 2.1, and the values of < N(s) > imputed 
from experiments [12], [13] and [14], and also observing 
that ^(s) is obtained from Eq. (12), we have estimated 
the values of /3 at various energies. The results are dis- 
played in Table 1 . We can see clearly that /3 increases as 
the collision energy also increases. /3 can be parameter- 
ized as Pis) = 77A8y/s + 0.4168 with x^/Ndf=1- 

Tabic 1 

/9(s) estimated values at various energies. The values of 
< N{s) > was imputed from Refs. [12], [13] and [14]. 



GeV 




< N{s) > 


I3{s) GeV 


52.6 


1.612 


11.55 


69.295 


200 


1.517 


21.4 


203.531 


546 


1.424 


27.5 


292.729 


900 


1.377 


35.6 


452.376 


2360 


1.286 




1061.13 


7000 


1.188 




2995.08 


14000 


1.130 




5912.68 



Now we proceed to compute the impact parameter de- 
pendence of inelasticity and infer some information on 
its energy dependency. 

3.3 Results and Discussion 

Based on the connection between multiplicity distribu- 
tion and inelasticity, we have update and applied the 
One String Model formalism deriving an expression, Eq. 
(22) , which allows us to study the impact parameter and 
energy dependence of inelasticity. Adopting the Block et 
al. QCD-inspired parametrization for Xpp{s, b) [6] and 
by using the estimated values of /3(s), Table 1, we have 
applied the Eq. (22) by fixing the collision energy and 
computed the inelasticity as a function of b. We show, 
in Fig. 7, the results from our analysis. Naturally, the 
inelasticity decreases as a function of impact parameter, 
b. The inelasticity behavior is essentially the same at en- 
ergies 546 and 900 GeV. At energies 52.6 and 200 GeV, 
we can sec appreciable difference just in the region of 
6 ~ (central collisions). It is interesting to note that, 
from the range of collision energy 50 ~ 200 GeV to that 
one 500 ~ 900 GeV, the inelasticity shows a difference 
about 60 percent in its values for b < 0.5 fm. We can also 
see that, at fixed b, the inelasticity K decreases as -^/s 
increases in the interval 52.6 - 900 GeV. The Eq. (22) 
depends on the eikonal and /3 parameter. With the 
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Fig. 7. Impact parameter dependence of the inelEisticity by 
using the formula obtained in this work, Eq. (22). 

eikonal as determined phenomenologically in [6] as in- 
put, where high energy cross sections grow with energy 
as a consequence of the increasing number of soft partons 
populating the colliding particles {pp/pp), it seems quite 
natural to expect that multiparton interactions leads to 
larger multiplicities/inelasticities as consequence to the 
full development of the gluonic structure. However, look- 
ing the same impact parameter dependence of inelastic- 
ity functions at 546 and 900 GeV, Fig. 7, we would be 
tempted to conclude that we are observing saturation 
effects due gluon recombination in the inelasticity, but 
there is no evidence of saturation in this range of energy 
(52.6 - 900 GeV). We have also applied our approach at 
energies y^=2.36 and 7 TeV (LHC), as shown in Fig. 
7, and the results suggest that the inelasticity is an in- 
creasing function of energy for the interval 2.36 - 7 TeV. 
As mentioned before, our main purpose is study fea- 
tures of multiplicity distributions deriving information 
on inelasticity and our analysis is based on the model in 
which is assumed that in parton-parton collision there 
is formation of a string. Thus, despite some simplifica- 
tions made in the One String Model, the results seems 
to be consistent with the multiplicity distributions data 
in a wide interval of energy (52.6 - 900 GeV). Hence, 
the computed inelasticities, in this range of energy, are 
reliable results. In counterpart, the results at ■^5=2. 36 
and 7 TeV are inconclusive in the context our analysis, 
because the One String Model probably underestimates 
the high multiplicities events due to the lack of the mul- 
ticomponent structure in its formulation. In fact, recent 
results reported by CMS Collaboration pointed out the 
importance of a multicomponent structure in hadron- 
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hadron inelastic interactions, in agreement with previ- 
ous experimental results (for details see [15]). Inelastic- 
ity also has been studied recently in Ref. [19], where, in 
the context of both Wdowczyk and Wolfendale model 
and UHECR data analysis, it was found that the inelas- 
ticity decreases in very high energy interactions, and, 
in the same work and by using the modified Feynman 
scaling formula, the inelasticity is an increasing function 
of the energy. It reflects the subtlety of the theme. We 
note that at ISR Energies (30-60 GeV) , where the lead- 
ing particle spectrum could be measured, the inelastic- 
ity is defined to be about 0.5. This value can be iden- 
tified with pp collision taking place at 6 ~ 0.6 fm, Fig. 
7. The One String Model has been used to study the in- 
fluence on $ considering possible values of A parameter 
at -^5=546 GeV, Fig. 6. In addition, wc have also com- 
puted $ at energies 52.6, 200 and 900 GeV using the 
different values of A (0.258, 0.198 and 0.11) and the re- 
sults, in each energy, are similar with that obtained in 
Fig. 6. Finally, we emphasize that the curves in Figs. 1-4 
has not been fit to data, except for the values of A and 
k (fixed) no experimental information about multiplic- 
ity distribution has gone into the calculation. Hence, the 
energy evolution of the multiplicity distributions, from 
ISR to Collider (30 ~ 900 GeV), is correctly reproduced 
by changing only the Overlap Function, Eq.(2), with- 
out changing the underlying elementary interaction, in 
agreement with what could be expected from QCD. 

4 Concluding Remetrks 

Being the impact parameter b an essential variable in a 
geometrical description of hadronic collisions, we have 
investigated the b dependence of inelasticity and also in- 
ferred some information on its energetic behavior. By 
using a geometrical model we have derived an expres- 
sion for K based on the hypothesis of connection be- 
tween multiplicity distribution and inelasticity. We have 
adopted the Block et at. model in our analysis, where the 
eikonal functions Xqq{s.h) and Xqg{s^b) are needed to 
describe the lower energy forward data, while Xgg{s, b) 
contribution dominates at high energy and determines 
the asymptotic behavior of cross sections. We believe 
that the same impact parameter dependence of the in- 
elasticity at energies 546 and 900 GeV can have impor- 
tant implications for the underlying gluon-gluon dynam- 
ics. In fact, we are testing our formalism by using the 
QCD Eikonal Model in which the gluon may develop a 
dynamical mass [20], [21]. At energies -^5=2. 36 and 7 
TeV the One String Model parametrization can not be 
tested, hence the results do not allow for any conclusion. 
Finally, the results suggest that there are relationships 
between the inelasticity and the eikonal function. 
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